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I address some recently raised issues regarding the time-parametrization dependence in stochastic 
descriptions of eternal inflation. To clarify the role of the choice of the time gauge, I show examples 
of gauge-dependent as well as gauge-independent statements about physical observables in eternally 
inflating spacetimes. In particular, the relative abundance of thermalized and inflating regions is 
highly gauge- dependent. The unbounded growth of the 3- volume of the inflating regions is found in 
certain time gauges, such as the proper time or the scale factor gauge. Yet in the same spacetimes 
there exist time foliations with a finite and monotonically decreasing 3-volume, which I demonstrate 
by an explicit construction. I also show that there exists no "correct" choice of the time gauge that 
would yield an unbiased stationary probability distribution for observables in thermalized regions. 



I. INTRODUCTION 

In most inflationary scenarios, the exit from inflation 
and the subsequent reheating do not occur everywhere 
at once. Some regions of the universe finish inflation 
and thermalize, while other regions continue inflating. 
Generic models of inflation predict that such inflating 
regions will be found at arbitrarily late times. This phe- 
nomenon was called eternal inflation 1 0, 0, 0] • The in- 
flating regions expand much faster than the thermalized 
regions, so heuristically one expects that most of the vol- 
ume of the universe is dominated by regions that have 
been inflating for a very long time and thus retain no 
memory of the initial state. The planet Earth is equally 
likely to be in any of the (infinitely many) thermalized re- 
gions, hence our present-day observations are essentially 
independent of the conditions before the onset of infla- 
tion. In this way the picture of eternal inflation was found 
to alleviate the problem of initial conditions for inflation. 
Any initial state is permissible (agrees with our present 
observations) as long as there is a nonzero probability for 
eternal inflation to set in. 

The inflating (i.e. not yet thermalized) domain in an 
eternally inflating spacetime can be described as a ran- 
dom collection of horizon-sized regions with locally de 
Sitter geometry. All regions expand at different local 
Hubble rates H(t,x), where x is a comoving coordi- 
nate and t is the proper time measured along comov- 
ing geodesies x — const. A region with expansion rate 
H evolves during one Hubble time At ~ H^ 1 into e? 
new horizon-sized regions having slightly changed expan- 
sion rates; some of these new regions may accidentally 
thermalize. Thus the evolution of the spacetime can be 
viewed as a realization of a branching diffusion process. 
One can compute probability distributions for various 
physical observables, such as the scalar field <f> or the scale 
factor a, at a given time t. Calculations are performed by 
modeling the branching diffusion process using Langevin 



Since it was not implied that inflation is eternal also to the past, 
a more precise designation would be "future-eternal inflation." 



and Fokker-Planck (FP) equations (see £| for a review 
and further references). Below I shall summarily refer to 
this description as the FP formalism. 

A description through an equal-time probability dis- 
tribution necessarily suffers from a dependence on the 
choice of the time variable t and in some cases this gauge 
dependence is quite dramatic 0, IU, El- The source of 
the gauge dependence is a bias introduced into the se- 
lection of observers when a particular equal-time surface 
is chosen. For example, the proper-time gauge favors re- 
cently thermalized regions, as compared with the scale- 
factor gauge. Probability distributions can be defined in 
a gauge-invariant manner in some cases, for instance in 
a class of models where an observable assumes all pos- 
sible values within one causally connected thermalized 
patch 0,0. However, a straightforward gauge- invariant 
definition of probability distributions is still lacking in the 
important case where observables assume different values 
in causally disconnected patches (for some approaches, 
see ({J)- 

In this paper I shall not consider these gauge-invariant 
prescriptions but instead focus on the issues of gauge de- 
pendence in the usual FP formalism. This apparently 
unavoidable gauge dependence has caused some conster- 
nation jig . regarding the overall validity of the stan- 
dard picture of eternal inflation. According to a criticism 
put forward in Ref. inflating regions cannot be said 
to dominate the volume of the universe because the vol- 
ume is not a gauge-invariant characteristic. One might 
then ask whether the entire picture of eternal inflation is 
not a gauge artifact or a result of inconsistently applied 
statistics. For instance, it was found in Ref. that cer- 
tain volume-weighted observables behave quite similarly 
in models with and without eternal inflation, thus casting 
doubt on the usually claimed ability of eternal inflation 
to forget initial conditions. 

The purpose of this paper is to answer these criti- 
cisms and to describe the defining features of eternal in- 
flation more accurately. As I will discuss below, the men- 
tioned criticisms correctly point out the flaws in certain 
often made statements about eternal inflation. Neverthe- 
less, the standard picture of eternal inflation remains un- 
changed. The presence of eternal inflation is manifested 
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by the permanent increase of the number of causally dis- 
connected non-thermalized regions. On the other hand, 
answers to certain other questions, especially those in- 
volving the proper 3-volume, turn out to be heavily in- 
fluenced by the choice of equal-time surfaces. I shall at- 
tempt to clarify the influence of the gauge by analyzing 
various specific statements that might be made about 
observables in eternally inflating spacetimes. 
The following questions will be answered: 

• Whether the volume of inflating regions is larger 
than the volume of thermalized regions. (Depends 
on the choice of the time gauge.) 

• More generally, whether a volume-weighted value 
of an observable has a qualitatively different be- 
havior when eternal inflation occurs. (Depends on 
the choice of the time gauge.) 

• Whether it is possible to describe the manifesta- 
tions of eternal inflation in a gauge-invariant man- 
ner. In particular, whether the volume of the in- 
flating domain exhibits an unbounded growth with 
time regardless of the time gauge. (Yes, no.) 

• Whether the FP formalism admits a particular time 
gauge such that the resulting probability distribu- 
tion is unbiased. (No.) 

II. ETERNAL INFLATION IN A BOX 

The global geometry of eternally inflating spacetimes 
possesses certain counter-intuitive properties. To render 
the discussion more visual, let us turn to a drastically 
simplified toy model of eternal inflation that nevertheless 
exhibits all the qualitative features of interest. 

In the toy model which we call "inflation in a box," 
time elapses in discrete steps s = 0, 1, 2, ... and the space 
is reduced to two dimensions and further to the square 
domain < x, y < 1, where x,y are the comoving co- 
ordinates measured in Hubble units H . The entire 
initial Hubble-size domain is assumed to be inflating at 
time s — 0. To imitate inflation during one timestep, we 
subdivide the initial inflating square into N x N equal 
sub-squares of size TV -1 x TV -1 ; at the next step, each 
sub-square will again have the Hubble proper size. Then 
we randomly mark some of the smaller squares as "ther- 
malized" assuming that each Hubble-size inflating square 
continues inflation with a probability q (where < q < 1) 
and thermalizes with probability 1 — q, independently of 
all other squares. The selection of thermalized squares 
concludes the simulation for one timestep. At the next 
timestep, the same procedure of subdivision and ran- 
dom thcrmalization is applied to each Hubble-sized in- 
flating square, while the "thermalized" squares do not 
evolve any further (see Fig. ^) . This random process was 
called a "random Sierpihski carpet" or a "curdling" in the 

book H2. 




| inflating regions 
| thermalized regions 



Figure 1: First steps in the construction of a random Sier- 
piriski carpet with N — 5 and q — 5/6. 



At a step s we thus have a number of inflating squares 
with comoving size N~ a x N~ s and also a number of ther- 
malized squares whose comoving sizes are determined by 
the times of their thermalization. It is straightforward 
to show that a given comoving point (x, y) remains in 
the inflating regime with probability q s after s steps, and 
that the mean number of inflating squares at a step s is 
(N 2 q) s while the mean comoving 2- volume of the inflat- 
ing domain is q s . Since each inflating square has Hubble 
proper size, the number of inflating squares (N 2 q) is 
naturally interpreted as the proper (as opposed to co- 
moving) volume of the inflating domain at time s. It 
follows that the mean proper volume of inflating domain 
grows without bound if N 2 q > 1 and decreases to zero 
if N 2 q < 1. Thus the condition N 2 q > 1 allows future- 
eternal inflation to occur with a nonzero probability. (A 
rigorous derivation of these results is given in the theory 
of branching processes, see e.g. the book 0] , chapter 1.) 

This toy model serves as a crude analogy for the pro- 
cess of random nucleation of bubbles in a de Sitter space- 
time (see e.g. 0|). The "thermalized squares" represent 
the nucleated bubbles of true vacuum which expand and 
quickly reach the comoving Hubble size (computed at nu- 
cleation time). The subdivision of the inflating squares 
mimics an exponential expansion, while the parameter N 
plays the role of the expansion factor during one timestep. 
Since the thermalized squares are not further subdivided, 
the expansion of spacetime within the true vacuum bub- 
bles is not modeled. As a remedy we consider a variation 
of the box model where the thermalized squares remain 
thermalized but are subdivided into M x M parts at each 
timestep, where M < N is chosen to emulate the slower 
expansion of true vacuum regions. This variation of the 
box model characterized by the parameters q, N, M will 
be sufficiently general for the present analysis. 

Let us now make a connection of the box model with 
models of eternal inflation driven by a scalar field 4> with 
an effective potential V(tfi). The FP formalism requires to 
choose a time slicing parametrized by a time coordinate 
r. A possible choice is r = t where t is the proper time 
measured along comoving geodesies x = const. Then one 
considers two probability distributions defined at a fixed 
time r: first, the distribution P(</>, r) of values along a 
given (randomly chosen) comoving trajectory, and sec- 
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ond, the total proper volume Pv(<fi, t) of all regions with 
the value <f> at time t. (Some issues regarding the inter- 
pretation of the distributions P and Py are clarified in 
Appendix [S]) If the time parameter r is related to the 
proper time t by a function of <j) only, 

dr = T{(j))dt, T((p) > 0, (1) 

then one can compute the probability distributions P and 
Pv by solving the FP equations, 2 

d T P = d 2 (D(4>)P) - (v(cf>)P) , (2) 

d T P v = dl {D{ct>)Pv) - 9 {v{<t>)Pv) + 3h(cf>)P v , (3) 

where Did), v (<t>), and h(<j)) are appropriate kinetic coef- 
ficients [la,ll3 that depend on the potential V((j>) and on 
the function T(<fi). The precise form of these coefficients 
will not be used in this paper. 

We can describe "inflation in a box" by analogous 
equations if we assume that the variable cj> takes on a 
discrete set of values, e.g. <j) = standing for infla- 
tion and 0=1 for thermalization. The distributions 
P(<p,r) and Pv{<t>, t) are reduced to Pj{s) = P{j,s) and 
Pvj(s) = Pv(j, s), where j = 0, 1 and s is the timestep. 
The equations describing these distributions are 

P (s + l)=qP (s), (4) 

P 1 (s + l)=P 1 (s) + (l-q)P (s), (5) 

Pvo(s + 1) = N 2 qP V0 (s), (6) 

P vl (s + 1) = M 2 P vl (s) +N 2 (1- q) P vo (s), (7) 

where q, N, M are the parameters of the box model. With 
the initial conditions P (0) = Py (0) = 1, Pl(0) = 
-FVi(O) = 0, the solutions are 

q s ,Pvo(s) = (N 2 q) S , P 1 («) = l-g« 
1-9,.,, 1S 1-A S . M 2 



With the same initial conditions as above, the solution is 



Po(s) 

Pvi(s) =- 



(N 2 q y 



A = 



(8) 
(9) 



q ~ 1-A' " N 2 q' 

It is clear that the proper volume of the inflating domain 
Pvo(s) grows without bound if N 2 q > 1. This is the 
usual sign of eternal inflation. 

It will be more convenient to analyze the box model 
in continuous time t rather than in the discrete time s. 
Denoting by a the thermalization rate and by Hq, Pi 
the expansion rates of the inflating and the thermalized 
domains respectively, we can write the following equa- 
tions describing the distributions of the comoving and 
the proper volume, 



dt 



P = -aP Q , 



dt 



Pi = aPn 



dt 
d 

dt 



Pvn = 2HnP V n - OR 



vo, 



P V1 = 2H 1 P vl +aP ] 



vo- 



(10) 

(11) 
(12) 



2 Here and below we use the Ito factor ordering in the diffusion 
terms 1 1 -Til - 



P (t)=e- at ,P V0 (t) = e 
P 1 (t) = l-e- at ,P V i(t) 



_ (2H -a)t 



3 (2ff -a)t _ p 2Hit 



(13) 
-a. (14) 



2(P -Pi) 

The correspondence with the discrete-time model is ob- 
tained by setting q = e~ a and Hq = In TV, Hi = InM. 
Eternal inflation is possible if a < 2Hq. 

An example of a different choice of the time variable is 
the "e- folding time" r = In a, where a is the scale factor. 
To derive the analogs of Eqs. l|lUD - (|12fl for the volume 
distributions at constant r, we note that an infinitesi- 
mal interval dr corresponds to the proper time interval 
dt = H ~ 1 dr within the inflating domains, and thus a 
fraction adt = aH^ 1 dr of the total inflating volume will 
thermalize during the interval dr. This applies to the co- 
moving as well as to the proper volume since the inflating 
domain is homogeneous. Therefore the increase <Ap in 
the comoving volume of the thermalized domain during 
the e-folding time interval dr will be aH^ 1 P dr. In this 
way we arrive to the equations 

■^-P VQ = 2P V0 - -^P yo , (16) 
dr Hq 

—P vl = 2P V1 + -^-Pvo- (17) 

Note that these equations cannot be obtained from 
Eqs. (|10|I - H12|I by a simple change of the time variable. 
The solutions of Eqs. ifTS Jl -lfTT jl are 



P (r) = e- aT ' Ha , P vo (t) = e (2-<*/H«)^ 



Pi{T) = l-e- aT ' H \P vl (T)=e 



- p2r _ (2-at/flo) 



(18) 

\ (19) 



The volume of the inflating domain grows proportion- 
ally to both the time gauges because the in- 
flationary expansion rate Ho is everywhere the same. A 
two-dimensional domain whose proper 2-volume grows 
slower than a 2 , namely as a 7 with 7 < 2, can be visual- 
ized as an expanding lacunary fractal with dimension 7, 
in the following sense. A lacunary fractal set S having a 
fractal dimension 7 < 2 contains infinitely many "holes" 
of diminishing size, and the 2-volume of S vanishes. To 
obtain a domain with a nonvanishing volume, one can 
coarse-grain the set Sona fixed scale which is the Hub- 
ble scale in the context of inflation. We may denote the 
result of the coarse-graining by S; the set S is interpreted 
as the inflating domain and the holes in S are regions that 
will eventually thermalize. At a given time, holes in S 
having sizes below the Hubble scale are invisible in the 
coarse-grained set S. As the set S expands (while the 
coarse-graining scale is kept fixed), it appears that new 
holes are constantly created within the domain S, mak- 
ing its volume grow slower than a 2 . The resulting fractal 
structure of the inflating domain was first explored in 
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Ref. |18| where it was shown that a domain expanding as 
oc a 1 is characterized by the fractal dimension 7. Thus 
the fractal dimension of the inflating domain in the box 
model is 7 = 2 — ctH^ 1 . A more rigorous definition of the 
fractal dimension of the inflating domain, together with 
a proof of its gauge independence, was given in Ref. [T^ . 



III. VOLUME OF INFLATING VS. 
THERM ALIZED DOMAINS 

Let us now consider the question of whether the inflat- 
ing domain dominates the volume of the universe at late 
times. In the box model with continuous time, the ra- 
tio of the total thermalized volume to the total inflating 
volume at a proper time t is found from Eqs. (|13[l - i|14[) . 



therm 



Pvijt) 
iVo(t) 



1 



-fit 



li = 2{H -H 1 ) 



(20) 

At late times t -> 00, this ratio tends to infinity if fi < 
and to a finite nonzero limit a/fi if fi > 0. The value 
of the latter limit may be any real number depending on 
the values of the parameters H ,Hi,a. With the typical 
assumptions a <C 2Hq and Hi < Hq, we find a finite and 
small value of the limit, a/fi <C 1, indicating that most 
of the volume is in the inflating domain. 

However, if we compute the volume ratio at a fixed 
scale factor rather than at a fixed step, we obtain a dif- 
ferent result, 



^thcrm 



V infi 



Pvi(r) 



Pvo(t) 



1. 



(21) 



This ratio tends to infinity as r — > 00, indicating that 
most of the volume is contained in the thermalized do- 
main, in a direct opposition to the typical result obtained 
in the proper time gauge. Therefore the statement that 
the volume is dominated by the inflating domain during 
eternal inflation does not hold in all gauges. Also, the 
ratio of the inflating to the thermalized volume does not 
necessarily characterize the presence or the absence of 
eternal inflation. 



the field <fi serves only to indicate the inflating or the 
thermalized state of a region, a volume-weighted average 
of 4> or, more generally, of a function A(<f)) at a time t is 
completely determined by the distributions Pvo(t) and 
Pvi(t), i.e. 



(A(<f>)) v = 



A(0)P VQ + A(1)P V1 _ A(0) + A(l) 



Pvo 



P\ 



VO 



Pvi 



1 



Pvi 
Pvo 



(23) 

As we have seen, the ratio Pvi /Pvo is highly gauge- 
dependent and so are the volume-weighted averages of 
<j) or of any functions of (/). 

Turning now to the scale factor a, we first note that in 
the scale factor gauge the entire equal-time hypersurface 
has the same value of a and it is useless to perform a 
volume-weighted averaging. Therefore we consider the 
proper time gauge, where the volume-weighted average 
of the scale factor differs from its statistical mean. To 
describe the distribution of the scale factor observed at a 
proper time t at a given comoving point in the box model, 
it suffices to compute the probability density Pth{s) for 
the given point to thermalize within a time interval [s, s + 
ds], 



Pth(s)ds 



s ads. 



(24) 



It is clear that pth(s) is the also the probability density 
for the scale factor to assume the value a = e H os+H 1 {t-s) 
at time t > s at a randomly chosen comoving point. It is 
more convenient to perform computations with the num- 
ber of e-foldings, In a, for which we have the statistical 
mean 

(In a) = (sH Q + (t - s) Hi) = H x t + (H - Hi) (s) (25) 

and the volume-weighted mean 

(lna) v = H 1 t+(H -H 1 )(s} v , (26) 

so it remains to compute the means (s) and {s) v using 
the distribution pth(s). Calculations yield 



Pth(s)sds 



1 - e" 



te~ 



(27) 



IV. VOLUME- WEIGHTED OBSERVABLES 

We now consider volume-weighted observables in an 
eternally inflating spacetime. A volume-weighted observ- 
able is defined using the volume V(A) occupied by regions 
that have a certain value A of the observable at a time t: 



E A AV(A) 
EaV(A) 



(22) 



Volume- weighting is always performed along a particular 
equal-time surface. 

In the box model, the possible candidates for volume 
weighting are the field and the scale factor a. Since 



and 

_ JoP th (s)e 2H " s+2H ^- s hd S _ pte* - e"' + 1 
JoPth(s)e 2H « s + 2H ^- s )ds ~ (e^-l)/x ' 

(28) 

where p is defined in Eq. (|20|l . Depending on the value 
of the parameter /i, the volume- weighted average for late 
times t may be dominated by the contribution of either 
inflating or thermalized regions. For instance, if fi < 
(domination by thermalized regions) then 



lim (s), 



1 



(29) 
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and thus for large t we have 

{lna) v = H 1 t--(H -H 1 )^H 1 t, (30) 
A* 

which is almost the same as the statistical mean, 

(hxa)=H 1 t + -(H -H 1 ). (31) 
a 

On the other hand, if fj, > (domination by inflating 
regions) then (s) v f» t, i.e. almost all volume is filled by 
very recently thermalized regions, and so 

(lna) v . « ff^ + (H - Hi)t = H t. (32) 

Note that the possibility of eternal inflation is compat- 
ible with both cases fj, > and /i < 0. We conclude 
that volume- weighted observables such as (lna)y gener- 
ally furnish neither any information specific to eternal 
inflation nor any gauge-invariant information. 

V. GROWTH OF VOLUME OF INFLATING 
DOMAINS 

The hallmark of eternal inflation is the unbounded in- 
crease in the total number of independent inflating re- 
gions. It is also often stated that the total proper 3- 
volume of the inflating domain grows without bound at 
late times in an eternally inflating spacetime; this is cer- 
tainly true for the 3-volume computed along hypersur- 
faces of equal proper time or of equal scale factor. One 
of the criticisms expressed in Refs. 0, |n| was that the 
3-volume is a gauge-dependent characteristic and thus 
cannot be used as an indication of the presence or the 
absence of eternal inflation. In an earlier attempt to re- 
solve this issue, it was proven in Ref. 0] (Sec. Ill A) that 
an eternally inflating spacetime manifests an unbounded 
volume growth in all time gauges of the form Q . How- 
ever, one may imagine time gauges not of the form ffl. 
Equal-time probability distributions in such gauges can- 
not be described by the FP equations @-@ where the 
kinetic coefficients depend only on the field (p. 

In fact, the statement about the growing volume of 
the inflating domain is gauge-dependent. If one allows 
completely arbitrary time slicings, the total volume of 
the inflating domain during eternal inflation might either 
grow or decrease depending on the time slicing. I illus- 
trate this dramatic dependence on the choice of gauge in 
Appendix FBI where I present an explicitly constructed fo- 
liation of a de Sitter spacetime by spacelike hypersurfaces 
whose proper 3-volume monotonically decreases to zero 
at late times. These hypersurfaces are artificially cho- 
sen to be almost null almost everywhere, which makes 
their 3-volume arbitrarily small despite the expansion of 
the background spacetime. Since an inhomogeneous eter- 
nally inflating universe expands "slower" than a pure de 
Sitter universe, the same family of hypersurfaces will also 
have a decreasing 3-volume in inhomogeneous inflating 



universes. It follows that the 3-volume of an arbitrary 
family of equal-time hypersurfaces cannot be used as a 
criterion for the presence of eternal inflation. 



However, a weaker statement is sufficient: namely, 
eternal inflation is present if (and only if) there exists a 
choice of time slicing with an unbounded growth of the 3- 
volumc of inflating domains. In the remainder of this sec- 
tion we shall demonstrate the equivalence of this criterion 
and the usual definition of eternal inflation (the existence 
of inflating domains at arbitrary late times). Note that 
the existence of something "at an arbitrarily late time" is 
a gauge-independent notion due to the monotonicity of 
time parameters in all gauges. 



If the 3-volume of the inflating domain grows without 
bound in some time gauge, then obviously there exist 
inflating domains at arbitrarily late times. It remains 
to prove that if inflating domains exist arbitrarily late, 
then the inflating volume grows in some gauge. This 
growth can be most easily demonstrated in the e-folding 
gauge, r = In a, using the construction of "eternal comov- 
ing points" 01 • These points follow comoving geodesic 
worldlines that forever remain within the inflating do- 
main and never enter any thermalized regions. It was 
shown in Ref. 0] using general topological arguments 
that the presence of inflating domains at arbitrarily late 
times entails the existence of infinitely many such "eter- 
nal points." Since each eternal point is surrounded by 
at least a Hubble-sized inflating region at any time, the 
3-volume of the inflating domain at e-folding time r is 
greater than H~ 3 N, where N is the number of eternal 
points that are pairwise separated by a proper distance 
of at least H^ 1 at time r. To show that the volume is 
greater than a preassigned bound M , we thus need to 
find at least No = MH 3 widely separated eternal points 
at a sufficiently late time. Since the eternal points are 
infinitely numerous, we can always choose a subset S 
containing No such points and then determine the mini- 
mum comoving distance Sx between any two points from 
the subset S. Within a hypersurface of constant scale 
factor a — e T , the proper Hubble distance H^ 1 cor- 
responds to the comoving distance H~ 1 e~ T . For suffi- 
ciently late e-folding times r > tq = — In(HSx) the scale 
factor a = e T is large enough so that the proper distance 
between any two points from S is larger than H~ x . Thus 
for t > To there exist at least Nq independent inflating 
regions, each having a volume at least H~ 3 . So we have 
shown that the volume of the inflating domain for t > To 
is greater than NoH~ 3 = M. Since M is an arbitrary 
bound, we conclude that the 3-volume computed in the 
e-folding gauge grows without bound. It follows that an 
unbounded growth of the 3-volume is found also in all 
gauges of the form including the proper time gauge. 
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Figure 2: Illustrative inflationary potential with a flat self- 
reproduction regime fa < < fa and deterministic regimes 



< 6 < fa and 02 < 



IS THERE A "CORRECT" TIME 
PARAMETER? 



VI. 



As we have seen in Sec. IIIII the division of 3- volume 
between the inflating and the thermalized domains is 
strongly gauge-dependent. For instance, the proper-time 
distribution is biased towards regions with faster expan- 
sion rates, while the scale factor gauge favors regions 
with larger comoving volumes. The probability distri- 
bution P(0, t) necessarily depends on the choice of the 
time parameter r. It has been sometimes proposed that a 
certain choice of the parameter r will yield physically jus- 
tified, unbiased probability distributions. In this section 
I shall demonstrate that no choice of the time parameter 
guarantees unbiased results, at least in a certain class of 
inflationary models. 

In the relevant class of models, the self-reproduction 
regime is symmetric with respect to the observable pa- 
rameters; a specific model of this kind was considered in 
Ref. 0- For the present purposes it is sufficient to ana- 
lyze a toy model of scalar-field inflation with an effective 
potential shown in Fig. [21 Our analysis shall partially 
follow Ref. where a similar potential was introduced. 
The potential is flat in the range 0i < < 2 where 
the evolution is fluctuation-dominated, while the evolu- 
tion of regions with > 02 or < ^i is completely 
deterministic (fluctuation- free). It is assumed that the 
fluctuation-dominated range fa < < fa is sufficiently 
wide to cause an eternal self-reproduction of inflating re- 
gions. There are two thermalization points, <f> — 0* 

(2) 

and = 0* , which may be associated to different types 
of true vacuum and thus to different observed values of 
cosmological parameters. The question is to compare the 
volumes Vi l2 of regions thermalized into these two vacua. 
Since there is an infinite volume thermalized into either 
vacuum, one hopes to obtain a sensible answer for the 
volume ratio V1/V2. 

In the present case, due to the symmetry of the po- 
tential it is natural to assume that Hubble-sized regions 
exiting the eternal self-reproduction regime at <f> = fa 
and at = 02 are equally abundant. Since the evolu- 
tion within the ranges 0^ < <f> < fa and 02 < 4> < ^l 2 ' 
is deterministic, the regions exiting the self-reproduction 
regime at = fa or = fa will experience a fixed num- 



ber of e-foldings which we may denote InZi and ln^ 
respectively. Therefore the volume of regions thermal- 



zed at = 



k<J) 



where j — 1,2, will be increased by the 



factors Z? and thus the volume ratio is 



Vi 

V 2 



z? 

zr 



(33) 



We shall now compare this result with that of the 
gauge-dependent FP approach. The time-dependent dis- 
tribution of the 3- volume Py (0, t) is a solution of the FP 
equation Q in a time gauge r. The late-time behavior 
of the volume distribution is 



P 



pf (<p)e^, 



(34) 



where 7 > is the largest eigenvalue of the corresponding 
stationary FP equation. (It was shown in Ref. Jljj that 
7 > in all gauges r admitted by the FP formalism.) 
Since the volume of regions with all values of 4> grows 
without bound, one needs to introduce a cutoff to com- 
pare the total volumes Vi,2 thermalized at <fi = cf>* and 
(2) 

<p = <p* ■ The equal-time cutoff, i.e. the procedure where 
one counts only the volume of the regions thermalized 
before a fixed time r = r max , is implemented straight- 
forwardly in the FP formalism and in the present case 
yields 



■dry 
6] 



^ = ^exp[-7(An-AT 2 )], (35) 



where Ar, , j — 1,2, are the time intervals elapsed in the 
time gauge r during the deterministic evolution of the 



field from <j)j to 



One can derive Eq. I|35l) from elementary consider- 
ations as follows. A Hubble region exiting the self- 
reproduction regime at <p — fa will thermalize after 
a fixed duration of time An and will accumulate the 

growth factor Z\ during this time. Similarly, regions exit- 

(2) 

ing at = 02 will reach the thermalization point = 0* 
after a time interval At2 and accumulate the growth fac- 
tor Z2- The "arrival times" Ar, and the growth factors 
Zj, j = 1, 2, can be computed using the particular shape 
of the potential V(0) but we shall not need their explicit 
forms. 

(7) 

It follows that the regions thermalized at = 0* 
at time r have previously exited the self-reproduction 
regime at time r — Ar, . According to Eq. I|34() , the vol- 
ume of regions exiting at = fa at time r — Atj is 



p(° ) (0 J )exp hir-Arj)} 



(36) 



where P, 



v due to the symmetry 



r ,0 l) =pr(0 2 )^p 1 (o) 

of the potential in the self-reproduction regime. The total 
volume of regions thermalized at = 0l J ' up to a time 



( T max) — Py^ Zj 



exp[7(r- Arj)] 



oc Zf exp [7(r max - At,-)] , 



(37) 
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where we have omitted j-independent factors. Thus we 
obtain the volume ratio as the limit 

Vi Vi(w) z\ 

v" = hm ^Tu V = ^3 ex P -7(An - Ar 2 ] . 

V 2 T max ^oo l/ 2 (r max j 

(38) 

This coincides with Eq. QHfijl. 

As already noted in Ref. [6j, the reason for the gauge 
dependence of the volume ratio (|35() is the presence of the 
times Atj which generally leads to an appreciable bias. 
The arrival times Arj are generally different, An ^ Ar 2 , 
unless the inflaton potential V{4>) is completely symmet- 
ric. 3 Since 7 > in all gauges, it is impossible to choose 
a gauge r in which the results l|33|l and l|35(l would coin- 
cide for all potentials. We conclude that there exists no 
"correct" gauge r that would always guarantee unbiased 
results for the equal-time cutoff procedure. 
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Appendix A: INTERPRETATIONS OF THE 
FOKKER-PLANCK FORMALISM 

In this appendix I shall clarify the possible interpreta- 
tions of the distributions P(<p,r) and Pv(4>,t) in terms 
of an ensemble of imaginary comoving observers placed 
throughout the inflationary spacetime. 

The comoving volume distribution P(4>, r) can be 
found as a solution of the FP equation as well as 
via a Langevin equation. The Langevin and the Fokker- 
Planck equations are equivalent descriptions of the same 
stochastic process, and we shall focus on the FP descrip- 
tion. 

The distribution P((f), r) allows two interpretations. 
On the one hand, P((j), r) is the probability for the field 
to have a value <f> at time r along a single, randomly cho- 
sen comoving trajectory; on the other hand, P(4>, r) is 
the fraction of comoving volume where the field has the 
value 0, within an equal-r hypersurface. (The "comoving 
3- volume" of a subdomain S of an equal-r hypersurface 
is defined as the 3-volume of the subdomain So of the 
hypersurface r = 0, where So consists of points whose 
comoving worldlines pass through S at time r.) 



3 The arrival times would be equal if it were possible to use the 
hypersurfaces of equal </> as the time slicing. However, a hypersur- 
face (f> = const is in general not spacelike because the evolution 
of <j> is not monotonic in the self-reproduction regime. 



The equivalence of these interpretations can be seen 
from the following argument. Consider a large but finite 
domain of volume Vq within the hypersurface r = and 
a very dense grid of comoving trajectories starting out 
from this domain. The total number N of these trajec- 
tories may be chosen as finite but very large, and the 
initial volume Vq can be subdivided into N tiny regions 
of comoving volume Vq/N situated around each of the 
N comoving lines. If the initial volume Vq is sufficiently 
large so that all the statistically possible histories are suf- 
ficiently well sampled, then the number of lines that reach 
a value at a time r is very close to NP(<fi, r). Thus the 
comoving volume of the subdomain of the equal-r hyper- 
surface having the value <j> is VqP{4>, r). Therefore in the 
statistical limit N — > 00 and Vq — * 00, the fraction of the 
comoving volume containing the field value (f> at the time 
r is equal to P(<f>, r). 

Conversely, if P((j>, r) is the statistical distribution of 
the comoving volume along an equal-r hypersurface, then 
the comoving volume of the subdomain with the value <p 
will be VoP((f>, r). We denote this subdomain by S((j), r) 
and consider a comoving trajectory randomly chosen 
from the total comoving volume Vq. The probability of 
this trajectory passing through the subdomain S(<f>, r) is 
equal to V P{(j},T)/V = P(</»,r). Therefore P(4>,t) is 
the probability of the field having the value <f> at time r 
along a randomly chosen trajectory. 

We note that there do exist correlations between the 
values of <f> at nearby comoving trajectories. However, 
these correlations are not reflected by the distribution 
P{4>, r) since it specifies only the total volume, but not 
the location, of regions with the field value <f>. If the initial 
comoving volume Vq is sufficiently large (many Hubble 
volumes), the distribution of the field (f> at time r along 
a randomly chosen comoving trajectory will be P(<fi, r) 
despite the correlations between nearby trajectories. 

Turning now to the distribution of the proper volume 
Pv{4>,t), we note that there is now only one interpreta- 
tion, namely Py (<fi, r) is the total proper 3-volume of the 
subdomain with the field value <p within an equal-r hyper- 
surface. (As before, we consider only the comoving future 
of a large but finite initial volume Vq at time r = 0.) It 
appears to be impossible to interpret Pv{4>i T ) directly in 
terms of observations made on a single randomly chosen 
comoving trajectory. The distribution Py(4>, t) describes 
a branching diffusion process with a possibly unlimited 
growth, and its interpretation necessarily has to involve 
a potentially very large ensemble of trajectories. 

A connection can be formulated between the distribu- 
tion Pv(4>,t) and the joint comoving distribution of the 
field and the scale factor P((f>, a, r). The comoving dis- 
tribution P(4>, a, r) is interpreted as the fraction of the 
comoving volume where the field has the value <f> and the 
scale factor has the value a, within an equal-r hypersur- 
face. Alternatively, P{<j>, a, r) is the probability of the 
field having the value 4> and the scale factor having the 
value a at time r along a single, randomly chosen comov- 
ing trajectory. The distribution P{<p, a, r) is the solution 
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of the FP equation 

8 T P = L^P - d a {h{cp)aP) , (Al) 

where 

L^P = d\ (D(<j>)P) - 9 (v(<j>)P) (A2) 

is the differential operator entering Eq. (|2J) and h(<fi) is 
the "expansion rate" function defined by 

^ = hifia. (A3) 

For example, in the proper-time gauge h(<f>) = H(<j>) while 
in the e- folding gauge h = 1. 

To obtain the distribution Pv(4>,t), one can use the 
procedure of "volume- weighting," 

/>OC 

P v {4>,t)= I daa 3 P(ct),a,T), (A4) 
Jo 

which takes into account the expansion factor a 3 along 
the trajectory. It is easy to verify that the resulting (un- 
normalized) distribution Pv{<p,r) satisfies the FP equa- 
tion This establishes a link between the volume in- 
terpretation of the distribution Pv(4>,r) and properties 
of a single comoving trajectory embodied by the joint 
probability P{4>, a, t). 

Appendix B: COUNTEREXAMPLES TO THE 
GROWTH OF VOLUME 

In this section I shall present two examples of foliations 
described by spacelike hypersurfaces 

t = t(r,x), -oo < t < +oo, (Bl) 

such that the total proper 3-volume of each equal-r hy- 
persurface is finite and monotonically decreases with r, 
asymptotically approaching zero at late times. The two 
examples will be foliations of the Minkowski spacetime 
and of the future part of a de Sitter spacetime, the latter 
obviously interpreted as a homogeneous, eternally inflat- 
ing universe. These examples unequivocally demonstrate 
that an unbounded growth of the 3-volume of equal-time 
hypersurfaces during eternal inflation is not observed in 
all gauges. 

I shall first describe the idea behind the construction 
of the foliations and later give technical details. For sim- 
plicity, let us begin by considering a flat 1+1-dimensional 
spacetime with standard coordinates (t, x) and a saw- 
shaped hypersurface 

t = S a (x) = - cos" 1 fcos ^) , (B2) 

where the inverse cosine function is by definition such 
that < Sq(x) < L. The hypersurface t = Sq(x) is 
piecewise null and has kinks at points x — nL, n S Z 







t 










L 














X 


-2L 


-L 


L 


2L 


3L 



Figure 3: The piecewise-null surface So(x) having zero proper 
length. 

(see Fig. It is clear that the proper length (i.e. the 
"I-volume") of this hypersurface is equal to zero. We 
can then choose a smooth spacelike hypersurface situated 
sufficiently close to t — So(x), so that the proper length of 
the portion of that hypersurface corresponding to |i| < L 
is smaller than a chosen bound e. It is also possible to find 
another hypersurface for L < \x\ < 2L approaching the 
null surface So(x) even closer, so that the corresponding 
proper length is smaller than e/2, and so on. Thus there 
exists an unbounded spacelike hypersurface S £ (x) with a 
finite and arbitrarily small proper length C[S e ] < e. 

The Minkowski space can be foliated by copies of the 
surface t — Sq(x) shifted in the t direction. Instead of 
this "null foliation" we consider a family of spacelike hy- 
persurfaces labeled by a time parameter r, 

t = t (t,x) = T + S e ( T )(x), (B3) 

where e(r) is some function. The family (|J33|) will be an 
admissible foliation of the 1+1-dimensional Minkowski 
spacetime as long as dt/dr > and each hypersurface is 
spacelike, 

{d?-dx*)\ t=t[Tx)T=coQst <i). (B4) 

These conditions can be satisfied by a suitable choice of 
the function e(r), and moreover it is possible to have 
s(t) monotonically decrease to zero as r — > oo. The 
resulting family of hypersurfaces is the desired foliation 
of the Minkowski space with a finite and monotonically 
decreasing proper volume. 

This construction can be generalized to the 3+1- 
dimensional Minkowski spacetime by considering a family 
of spherically symmetric hypersurfaces 

t = t(T,r)=T + S e{T) (r), r=|x|. (B5) 

As will be demonstrated below, the function e(r) can be 
chosen to achieve a finite proper volume V(t) of the entire 
equal-r hypersurface and moreover to have V(r) — > as 
r — > oo. 

A slightly modified but basically similar construction 
yields an analogous family of spacelike hypersurfaces fo- 
liating a de Sitter spacetime. For simplicity we consider 
a 1 + 1-dimensional de Sitter spacetime with flat spatial 
sections. Then we can use conformal coordinates (rj,x), 



Figure 4: A null foliation of the de Sitter spacetime in con- 
formal coordinates (ji,x). Thicker lines are the consecutively 
smaller copies of the saw-shaped surface Sa(x), while the thin- 
ner lines in between show the continuous deformation proce- 
dure ("growing ridges"). The dashed line corresponds to an 
infinite future (77 — > 0). 



so that the de Sitter line element is conformally flat, 



ds 2 = dt 2 - e- 2Ht dx 2 = —5-; - (dr, 2 - dx 2 ) , (B6) 

H z r) 2 x ' 



-HI. 



H 



-00 < i] < 0. 



(B7) 



All lines 77 = 770 ± x are null geodesies and thus the saw- 
shaped hypersurface 77 = Sq(x) is piecewise null. How- 
ever, we cannot use the hypersurface 77 = Sq(x) to foli- 
ate the spacetime as in the Minkowski case because the 
boundary 77 = prevents us from shifting Sq (x) arbitrar- 
ily in 77. Instead it is possible to continuously deform the 
"saw" and to decrease the size of its ridges while grow- 
ing new ridges as time passes (see Fig. 0J. The resulting 
family of null hypersurfaces Sq(t,x) can be labeled by 
a conformal time parameter r, where r < 0, r — > as 
77 — s- 0, and it is clear from the figure that Sq(t, x) is a fo- 
liation of the entire future part of the de Sitter spacetime. 
The piecewise- null hypersurfaces So (t, x) are then used 
as a "skeleton" for a family of non-intersecting, smooth 
spacelike hypersurfaces 77 = H e (r, x) with a finite proper 
volume V[ff e ] < e. It is possible to choose the function 
s(t) such that e(r) — > as r — > 0. This procedure con- 
structs a foliation of the de Sitter spacetime with the 
total proper volume decreasing to zero at late times. 

In the remainder of this appendix I describe the re- 
quired hypersurfaces explicitly and verify their properties 
as claimed. 

The first step is to replace the angular piecewise-null 
hypersurface Sq(x) by a smooth hypersurface Sq(x;e) 
which is null everywhere except the "caps" at the cor- 
ners. To this end we use a "smooth corner" function 
C(x) which is infinitely differentiable, convex, and sat- 
isfies C(x) = \x\ for \x\ > 1 and \C'\ < 1 for |arj < 1. An 
example of a suitable function C(x) is plotted in Fig. [S] 



Figure 5: The "smooth corner" function C(x) defined by 
Eq. HHJ with A = 0.01, B = 0.5, Ci as 0.6613, C 2 « 3.323. 



and is defined for |x| < 1 by 

=1 



C(x) = l-d + C 2 







dy / dz exp 
'0 



A 

z 



D 



where A, B are positive constants and 

A 

z 



Ci 



f 1 f v 
C2 / dy dz exp 

Jo Jo 



D 



Co 



dz exp 



.4 

z 



D 



1-z 



1-z 

-1 



(B8) 



(B9) 
(B10) 



The function C(x) is smooth because by construction 
all derivatives of C[x) vanish at x = and at x = ±1 
[except C"(±l) = ±1], while at other x the function C(x) 
is analytic. Note also the inequality 



< C(x)-xC'(x) < 1 



(Bll) 



which follows from the convexity of C{x). 

A rescaled "smoothed corner" function C(x) will be 
used to replace small neighborhoods of the top and the 
bottom corners of the curve Sq(x). Namely, we define 
the "rescaled corner" functions 

C s (x) = C (xS- 1 ) S, C s (x) = L-C 5 {x), (B12) 

plotted in Fig. A spherically symmetric hypersurface 
t = Cg (|x|) is null except for a small range |x| < S, and its 
proper 3-volume is smaller than the volume of a sphere 
with radius S. We now replace the corners at \x\ — nL in 
the curve Sq(x) by smoothed corners Cg n (x) or Cs n (x), 
where 5 n decreases with n. For later convenience we set 



Sq — eL; S n — eLr 



71=1,2, 



(B13) 



where e < 1 is a free dimensionless parameter. To be 
definite, let us choose the smoothed-out version of Sq(x) 
as follows, 



S (x;e) = 
j C Sn (x-2nL), 



2n _I<M<2n+i; 
C 5n (x-(2n + l)L), 2n+ | < ^ < 2n+ |, 



(B14) 
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It then follows from Eq. (|B15|) that 



0<S e (r,r)-T<(l-h)L<L, (B20) 



dr 



dh 
dr 

Se(r,r) 



<l-h+L 



dh 



dr 



(B21) 
(B22) 



Therefore the hypersurface t — S £ (r, r) is everywhere 
spacelike. The proper 3-volume of this hypersurface is 



Figure 6: The rescaled "smooth corner" functions Cs{x) (cor- 
ner down) and Cg(x) (corner up) defined by Eq. 1B12II . 



where n = 0,1,2,... By construction, the function 
So(x;e) has the following properties that will be used 
below: 



0<S (x;e)<L, ^S Q (x;e) < 1; 



S (x;e) 



< sup 



dC s Jx) 



de 



<eL< L. 



(B15) 
(B16) 



The last property can be derived from Eq. 1|B11J| . 

The hypersurface t = >Su(|x ; e) is not everywhere null 
because of thin spacelike spherical shells in small neigh- 
borhoods of |x| = nL, n = 0,1,2,... The 3-volume of 
these shells can be estimated as 



A 00 



n=0 



< ( + 6i 2 ™ 2<J " ) < 30eL3 - ( B1? ) 



V [S e (r, r)\ = 4tt / dry/Wg(r), (B23) 



where is the determinant of the induced metric 
within the surface: 



ds 2 



(3) 



9(r) 



dS e 
dr 

dS £ 
dr 



dr 2 +r 2 dn 2 - 



(B24) 
(B25) 



We can now divide the hypersurface into the "caps," 
i.e. ranges of r within small neighborhoods of nL, n = 
0, 1,2, and "hills" where r is outside those neighbor- 
hoods. For r within the "hills" we have dSo/dr — ±1 and 
thus the following upper bound holds for those r, 



^dh dS 


dr 


dr 




dh 


+ 2 






dr 



Oh 



S < 2h + 2-L < Ah. 
L 



(B26) 



Therefore the total volume of all the "smoothing caps" 
can be made negligible by choosing a sufficiently small e. 

Now we are ready to describe the foliation of the 3+1- 
dimensional Minkowski spacetime. The spherically sym- 
metric hypersurface t — S e (r, r) is constructed by rescal- 
ing and shifting the smoothed function So(x;e) defined 
by Eq. (|B14() . A particular ansatz can be written as 

S e (r,r)=T+(l-h(r,r))So(r;£(r)), (B18) 

where h(r, r) and e(r) are suitable functions and < 
h(r,r) < 1. By construction, the surface S £ is every- 
where smooth (infinitely differentiable) and approaches 
the "skeleton" surface So as h — > 0, and we shall set 



h (t, r) = e exp 



2L 



e(r) = £ exp 



So the total volume of the "hills" is less than 
4tt^ r 2 y/4h(T,r)dr = 2048tt£L 3 exp (-^) ■ (B27) 

The total volume of the "caps" is constrained by 
Eq. (|B17|) . where e is now replaced by the function e(r) 
defined by Eq. (|B19() . Therefore an upper bound on the 
total volume is found as 



V [S e (T,r)\ < 30e(r)i 3 + 4vr / r 2 y/4h(T,r)dr 

Jo 



< 6464eL 3 exp 



(B28) 



(B19) 



As we intended, the total 3-volume of the hypersurface 
t — S e (r, r) is finite and approaches zero at late times 
(t — » oo). 

Finally we need to check that the hypersurfaces t = 
S e (r,r) are strictly monotonic in r, i.e. dS e /dr > 0. 
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Using Eqs. ljB16j> and l|Bl9j) . we find 



Figure 7: A typical piecewise-null hypersurface from the foli- 
ation in Fig. [I] 



(B29) 

(B30) 
(B31) 



dh 
Ih 
dS e 
dr 



dh„ ft 1 



1 



dh 

^rS Q + (i - h) 

OT 



dSo(r; e) de{r) 



de 



dr 



1 1 



Therefore the family of hypersurfaces S £ (r, r) is indeed a 
foliation of the entire Minkowski spacetime. This argu- 
ment completes the construction. 

Turning now to the future part of a de Sitter spacetime, 
say —H^ 1 < rj < 0, we start with the family of piecewise- 
null hypersurfaces depicted in Fig. A typical hyper- 
surface having smaller ridges of size / and larger ridges 
of size L is shown in Fig. We denote by S £ (x; I, L) a 
smoothed version of this hypersurface constructed using 
the "corner" functions (|B12|) ; a cumbersome explicit for- 
mula for S e (x\ l,L) can be written which we omit. One 
can easily verify the properties 



\S £ (x;l,L)\<L, 



dS F 



dx 



< 1. 



dS £ 



dS £ 



de 



< L. 



(B32) 

The desired family of smooth hypersurfaces rj = 
H £ (r,r) is built, as before, by shifting and rescaling the 
functions S s (x;l, L). An additional time dependence in 
the parameters I and L is needed to reproduce the be- 
havior indicated in Fig. 21 An explicit expression in con- 
formal coordinates (rj, r) can be written as 

n = H £ (t, r) = t + (1 - h(r, r)) S e(T) (r; l(r), L(r)), 

(B33) 

where we consider the range — H^ 1 < t < (recall that 
r is a conformal time parameter), the functions h(r,r) 
and e(r, r) are similar to those in the Minkowski case, 
and we define 

L(t)=2-^L, p(r) = - Llog 2 (-ifr)J , (B34) 

1(t) = 2- p( - t ^L+-HtL, (B35) 

where [^j denotes the integer part of x, i.e. the alge- 
braically largest integer not exceeding x. The defini- 
tions (|B34() - (IB35|) correspond to Fig. 0] and describe a 
gradual growth of ridges with size 1(t). By construction 
L{t) < L\Ht\; more precisely, during the time range 



_ h -i 2 -p+i < T < -H~ x 2-p we have L(t) = L2~p 
while the ridge size l(r) changes from to L2~ v . At 
t = —H~ 1 2~ p the ridge scale L is halved and I is reset 
to zero. This requires p > 0, i.e. r > —2H~ 1 , and the 
latter inequality holds since we consider only the range 
-iJ- 1 <r<0. 

The 3-volume of the hypersurface H £ (r, r) is found as 



V[H e (r,r)] = ^ I rdr 



y/1 - (dH £ /drf 

m 



(B36) 



Since the entire hypersurface is bounded in time, namely 
H £ (t, r) < |r, we may estimate its 3-volume from above 
using the maximum scale factor 8 \ Ht\ . An upper 
bound on the factor yl — {dH £ /dr) 2 is obtained simi- 
larly to Eq. (|B26() since the function S £ (r; I, L) has prop- 
erties <|B32I) similar to those of S (x;e), so 



1 



fdH £ 



\ dr 



< 4ft (r, r). 



(B37) 



We now choose e(r) and ft(r, r) as before except for the 
replacement r — > —H^ 1 In {—Ht) as appropriate for the 
conformal time, 



e(r) = e exp 



lnliJrl 



AHL 



ft(r, r) = e exp 



ln|lfr| r 
2HL ~ 2L 

(B38) 

Then the total volume of the hypersurface is estimated 
similarly to Eq. QB28I , 



V [H £ ] < 8 \Ht\ 



30e(r)i 3 



47r / r 2 ^Ah(r)dr 
/o 



< 51712eL J \Ht\* hl 



(B39) 



Choosing the free parameter L (the ridge size) such that 
HL < yj, i.e. having more than 12 ridges per Hubble 
length, we find 



MmV[H e (T,r)] = 



(B40) 



as required. 

It remains to verify that the hypersurfaces rj — H £ (r, r) 
are everywhere spacelike and monotonic in t. Note that 
the functions e(r) and ft(r, r) satisfy 



de 



e(r) 
AL \Ht\ 



< 



4L(r)' 



9ft 
y7 



ft(T, r) 
2L\Ht\ 



< 



2L(r) 
(B41) 

The spacelike character is proved by the estimate 



dH £ 



dr 



< 1 



which is derived from Eq. (|B32(I : 



d_ 

Or 



H £ (r,r) 



<(l-ft) 



S s (r,r) 



dh(r, r) 



di 



(B42) 



S e (r,r) 



<1 - ft- 



— L(t) < 1 
2L W " 



ft ft 

- = 1 < 1. 

2 2 

(B43) 
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The monotonicity, dH e /dr > 0, can be shown using 
Eq. i|ij32f> and the following inequalities, 



This argument completes the construction of the foliation 
for the de Sitter spacetimc. 



dh 
cV 



S e + {1- k) 



dS £ de{T) dS e dl{r) 



de dr 



dl dr 



(B44) 
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